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I I Abstract 

m ' Let 
> ; 

Au + k\{z)u = -^f{z) in X [0,1], (1) 
§■ m(x\0)=0, n'(x\l) = 0, (2) 

o : 

^ I where u = u{x^,z), := (xi, 2:2)1 f '■= := z, u' = (^(r) is the 

Oh' delta-function, n{z) is the refraction coefficient, which is assumed to be a real- 

r~| ■ valued integrable function, A; > is a fixed wavenumber. The solution to (l)-(2) is 

"j^ \ selected by the limiting absorption principle. 

^1 It is proved that if f{z) = 5{z — 1), then n(z) is uniquely determined by the 

^ ' data u{x^ ,1) known Vx^ G M^. Comments are made concerning the earlier study 
of a similar problem in the literature. 
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1 Introduction 

In [|l| the following inverse problem is studied: 



[A + kMz)]u = -^/(^), in X [0, 1], (l.i; 



^(a;^, 0) = 1) = 0, := (xi, X2), X3 := 2;, «':=—. (1.2) 



*Key words and phrases: inverse scattering, wave propogation, waveguides, ocean acoustics 
^Math subject classification: 35R30 



Here A; > is a fixed wavenumber, n{z) > is the refraction coefficient, which 
is assumed in [1] to be a continuous real-valued function satisfying the condition < 
n{z) < 1, the layer x [0,1] models shallow ocean, r := |a;^| = \/xi + x^, 6{r) is 
the delta-function, ^ = S{x^), f{z) G C^[0,1] is a function satisfying the following 
conditions m], p. 127: 

/(O) = /"(O) = /'(I) = 0, /'(O) ^ 0, /(I) ^ 0, f{z) > in (0, 1). (C) 

The solution to (1.1)-(1.2) in is required to satisfy some conditions ( p. 122, 
formulas (1.4), (1.8)-(1.10)) of the radiation conditions type. 

It is convenient to define the solution as u{x) = \imUi;{x), that is by the limiting 

absorption principle. We do not show the dependence on k in u{x) since k > is fixed 
throughout the paper. The function is the unique solution to problem (1.1)-(1.2) 

in which equation (1.1) is replaced by the equation with absorption: 

[A + k^n{z) - ie]u,{x) = -^f{z), in x [0, 1], e > 0. 

27rr 

One defines the differential operator corresponding to differential expression (1.1) and 
the boundary conditions (1.2) in L^(]R^ x [0,1]) as a selfadjoint operator (for example, 
as the Friedrichs extension of the symmetric operator with the domain consisting of 
H^iM."^ X [0, 1]) functions vanishing near infinity and satisfying conditions (1.2)), and 
then the function Us{x) is uniquely defined. By H"^ we mean the usual Sobolev space. 
One can prove that the limit of this function u{x) = limue{x) does exist globally in 

£j,0 

the weighted space L^(R^ x [0, 1], ^q^), a > 1, and locally in i/^(M^ x [0, 1]) outside a 
neighborhood of the set {r = 0, < 2; < 1}, provided Xj 7^ Vj, where Xj are defined 
in (1.7) below. This limit defines the unique solution to problem (1.1)-(1.2) satisfying 
the limiting absorption principle if Xj 7^ Vj. If f{z) = 6{z — 1), where 6{z — 1) is the 
delta- function, then an analytical formula for Ue(x) can be written: 

'^six) = ^i^j{z).fj—Ko{r^X] + ie), 
i=i ^ 

where KQ{r) is the modified Bessel function (the Macdonald function), and fj = '0^(1) are 
defined in (1.6) below, and i/^jiz) and are defined in formula (1.7) below. This formula 
can be checked by direct calculation and is obtained by the separation of variables. The 
known formula JF~^ A^+a^ ~ 2^-^o('2'") was used, and J-'u := u is the Fourier transform 
defined above formula (1.3). 

From the formula for u^{x), the known asymptotics Ko{r) = ^J^e~'^\V + 0(r^^)] for 
large values of r, the boundedness of |'0j(-2)| as j 00 and formula (1.8) below, one can 
see that the limit of u^[x) as e — does exist for any r > and z G [0, 1], if and only if 
Aj 7^ 0. If Xj = for some j = Jq, then the limiting absorption principle holds if and only 
if /jg = 0. If Aj 7^ \fj, then the limiting absorption principle holds and the solution to 
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problem (1.1)-(1.2) is well defined. If Xj = for some j = Jq, then we define the solution 
to problem (1.1)-(1.2) with f{z) = 6{z — 1) by the formula: 

11°° 1 
u{x) = ^j,{z)^j,{l)—\og{-) + ^j{z)i/jj{l)—Ko{rXj), r := \x^\. 

i=i.jyio 

This solution is unique in the class of functions of the form u{x) = Yl'jLi'^j{^^)'^j{^)y 
where AiUj-XjUj = -5{x'^) in M^, Aiw := w^^xi +Wx2X2, uj G L'^(R'^) if A| > 0; if A| < 
then Uj satisfies the radiation condition — i\Xj\uj) — > as r oo, uniformly in 

directions and if = then Uj = ^ log(^) + o(l) as r ^ oo. 

The inverse problem (IP) consists of finding n{z) given g{x^) := u{x^, 1) and assuming 
that f{z) = S{z-1) m (1.1). 

By the cylindrical symmetry one has g{x^) = g{r). 

It is claimed in p. 137] that the above inverse problem has not more than one 
solution, and a method for finding this solution is proposed. The arguments in [Q] are 
not satisfactory (see Remark 2.1 below, where some of the incorrect statements from [1], 
which invalidate the approach in [1], are pointed out). 

The aim of our paper is to prove that if f{z) = d{z — l), then n{z) can be uniquely and 
constructively determined from the data g{r) known for all r > 0. It is an open problem 
to find all such f{z) for which the IP has at most one solution. 

The method we use is developed in P (see also [^). Properties of the operator 
A + k'^n{z) in a layer were studied in p. In [§] an inverse problem for an inhomogeneous 
Schrodinger equation on the full axis was investigated. 

Let us outline our approach to IP. 

Take the Fourier transorm of (1.1)-(1.2) with respect to x^ and let 

y.-y{^z,X):=u:= f u{x\z)e''''<dx\ |C| := A, C e M^ 

and 

G{X) ■.= g{r). 

Then 

£v := v" — X'^v + q{z)v = —f(z), q{z) := k'^n{z), v = v{z,X), (1.3) 



t;(0,A) = t;'(l,A)=0, (1.4) 



t;(l,A) = G'(A). (1.5) 

IP: The inverse problem is: given G{X), for all X > and a fixed f{z) = d{z — 1), 
find q{z). 
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The solution to (1.3)-(1.4) is: 



^(^' ^) = E f^l' ^^'^ ^= M^A^)dz, (1-6) 

where ^^(2;) are the real- valued normalized eigenfunctions of the operator L :— —-^ — 
q{z): 

L^, = A>„ V',.(0)=V';.(1)=0, ||V',(^)|| = 1. (1.7) 

We can choose the eigenfunctions ipj{z) real- valued since the function q{z) = k'^n{z) is 
assumed real-valued. One can check that all the eigenvalues arc simple, that is, there is 
just one eigenfunction ipj corresponding to the eigenvalue A| (up to a constant factor, 
which for real- valued normalized eigenfunctions can be either 1 or —1). 
It is known (see e.g. [4. p. 71]) that 

A? = Aj - + 0{^)] as j - +00. (1.8) 
The data can be written as 

G(A) = E|^. (1-9) 

j = l j 

where fj arc defined in (1.6). The series (1.9) converges absolutely and uniformly on 
compact sets of the complex plane A outside the union of small discs centered at the 
points ±iXj. Thus, G(A) is a meromorphic function on the whole complex A-plane with 
simple poles at the points ±iXj. Its residue at A = iXj equals ^^2ix'^^ ■ 

If f{z) — 5{z — 1), then fj — ■0j(l) 7^ Vj = 1, 2, , (see section 2 for a proof of the 

inequality ■0j(l) 7^ Vj = 1, 2, ,) and the data (1.9) determine uniquely the set 

{Xl ^|(l)},=i,2,... (1.10) 

In section 2 we prove the basic result: 

Theorem 1.1. If f{z) — 5{z — l) then the data (1.5) determine q{z) e i^^(0, 1) uniquely. 

An algorithm for calculation oi q{z) from the data is described in section 2. 

Remcirk 1.2. The proof and the conclusion of Theorem 1.1 remain valid for other bound- 
ary conditions, for example, u'{x^,0) — u{x^, 1) with the data u{x^,0) known for all 
x^ e 
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2 Proofs: uniqueness theorem and inversion algo- 
rithm 



Proof of Theorem 1.1. The data (1.9) with f{z) = 5{z — 1), that is, with fj = 
determine uniquely {\^}j=i^2,... since ±iAj are the poles of the meromorphic function 
G{\) which is uniquely determined for all A e C by its values for all A > (in fact, by its 
values at any infinite sequence of A > which has a finite limit point on the real axis). 
The residues V'|(l) of ^(A) at A = iXj are also uniquely determined. 
Let us show that: 

i) V',(l)7^0 Vj = l,2,... 

ii) The set (1.10) determines q{z) G /^^(0, 1) uniquely. 

Let us prove i): 

If i'jiX) = then equation (1.7) and the Cauchy data '0i(l) = "^^(1) = imply that 
iljj{z) = which is impossible since || '>pj{z) ||= 1, where || u |p:= \u\^dx. 
Let us prove ii): 

It is sufficient to prove that the set (1.10) determines the norming constants 

:=|| ^^(2;) f 

and therefore the set 

{Aj, ttj}j=l,2,..., 

where the eigenvalues A^ are defined in (1.7), \E'j = ^{z, Xj), ipj{z) := -pqi^; 

_^"_s2^_g(2)^ = 0, ^(0,s) = 0, *'(0,s) = l, (2.1) 
and Xj are the zeros of the equation 

^'{l,s) = 0, s = X„j = l,2, (2.2) 

The function \E^'(1, s) is an entire function of u = s"^ of order i, so that (see 0): 

^'(l,s)=7j](^l-y , J = const. (2.3) 

From the Hadamard factorization theorem for entire functions of order < 1 formula (2.3) 
follows but the constant factor 7 remains undetermined. This factor is determined by the 
data {A|}vj because the main term of the asymptotics of function (2.3) for large positive 
s is cos(s), and the result in p. 243, (see Claim 1 below) implies that the constant 7 
in formula (2.3) can be computed explicitly: 
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where A° are the roots of the equation cos(s) = 0, A° = ^ j = 1^2, , and the 

infinite product in (2.3') converges because of (1.8). 

A simple derivation of (2.3'), independent of the result formulated in Claim 1 below, 
is based on the formula: 

1= hm 7-X = ^ll~6~- 

cosliw) 

For convenience of the reader let us formulate the result from 0], p. 243, which yields 
formula (2.3') as well: 

Claim 1: The function w{X) admits the representation 

riSin(A) h(\) 
w{X) = cos(A) - B — ^ + 

A A 

where B = const, h{X) = H{t) sm{Xt)dt, and H(t) E L'^iO, 1) if and only if 

°° A? - A2 



w 



where Aj = A° — y + y, Pj are some numbers satisfying the condition: J2'jLi l/^jP < cxd, 

Xj are the roots of the even function w{X) and A° = (j — i)7r, j = 1, 2, , are the positive 

roots o/cos(A). 
The equality 

where 7 is defined in (2.3'), is easy to prove: if w is the left-hand side and v the right-hand 
side of the above equality, then w and v are entire functions of A, the infinite products 

converge absolutely, ^-^^0)2 = ( 1 — 1 and taking the infinite product and using 

(2.3'), one concludes that ^ = 1, as claimed. 

In fact, one can establish formula (2.3") and prove that 7 in (2.3") is defined by (2.3') 
without assuming a priori that (2.3') holds and without using Claim 1. The following 
assumption suffices for the proof of (2.3"): 

i) A,2 = (A°)2 + 0(1), (A°)2 = 7r2(j - i)^ 

Indeed, if i) holds then both sides of (2.3") are entire functions with the same set of 
zeros and their ratio is a constant. This constant equals to 1 if there is a sequence of points 

at which this ratio converges to 1. Using the known formula: cos(A) = H^i (;^o)2 ; 

and the assumption i) one checks easily that the ratio of the left- and right-hand sides of 
(2.3") tends to 1 along the positive imaginary semiaxis. Thus, we have proved formulas 
(2.3)-(2.3') without reference to Claim 1. 
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The above claim is used with w{s) = in our paper. The fact that 

admits the representation required in the claim is checked by means of the formula for 
\l/'(l,s) in terms of the transformation operator: ^{z,s) = ^l^^l^fl -)- K{z, f) "'^^^^^ dt, 
and the properties of the kernel K{z,t) are studied in Thus, \E''(l,s) = cos(s) + 
■^(i.i)sm(s) _|_ ji t) ^^^^f^^ dt. This is the representation of \l/'(l,s) := w{s) used in 
Claim 1. 

Let us derive a formula for aj : = || |p. Denote \I/ := differentiate (2.1), with 
replaced by u, with respect to u and get: 

- z/^ - = ^. (2.4) 

Since q{z) is assumed real-valued, one may assume ■0 real-valued. Multiply (2.4) by 
\Ef and (2.1) by subtract and integrate over (0, 1) to get 



< a, := ^'^dz = (vl/;.vi/^. - 



-^.(imii), (2.5) 



where the boundary conditions ^^^(O) = ^l'j(l) = ^i(O) = were used. 
From (2.3) with = z/ one finds the numbers bj := ^^'(l): 



d 



A 
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Claim 2: The data "01(1) = — ^— ^j? where aj :=|| ^I/j(-2) o'l-'i equation (2.5) 
determine uniquely aj. 

Indeed, the numbers hj are the known numbers from formula (2.6). Denote by tj : = 
'?/'|(l) the quantities known from the data (1.10). Then it follows from (2.5) that = 
tjajb'j, so that 

a J = tjh]. (2.7) 

Claim 2 is proved. 

Thus, the data (1.10) determine aj =|| ^ j |p uniquely and analytically by the above 
formula, and consequently q{z) is uniquely determined by the following known theorem 
(see for example, [3]): 

The spectral function of the operator L determines q{z) uniquely. 

The spectral function p(A) of the operator L is defined by the formula (see [3, formula 
(10.5)]): 



A2<A ■ 



The Gelfand-Levitan algorithm [3] allows one to reconstruct analytically q{z) from 
the spectral function p(A) and therefore from the data (1.10), since, as we have proved 
already, these data determine the spectral function p(A) uniquely. 

Theorem 1.1 is proved. □ 
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Let us describe an algorithm for calculation of q{z) from the data g{x^): 

Step 1: Calculate G{\), the Fourier transform of g{x^). Given G{X), find its poles 

±iXj, and consequently the numbers A^; then find its residues, and consequently the 

numbers ^l)j{l)fj. 

Step 2: Calculate the function (2.3), and the constant 7 by formulas (2.3) and (2.3'). 
Calculate the numbers hj by formula (2.6) and aj by formula (2.7). Calculate the spectral 
function p(A) by formula (2.8). 

Step 3: Use the known Gel'fand-Levitan algorithm (see [3]-[5]) to calculate q{z) from 
p{\). 

This completes the description of the inversion algorithm for IP. 

Remark 2.1. There are inaccuracies in [Ij. We point out two of these, of which the 
first invalidates the approach in [Ij. 

In [1, p. 128, line 2] the an are not the same as an in formula [1, (3.3)]. If one uses 
an from formula [ 1, (3.3)], then one has to use in [1, p. 128, line 2] the coefficients 
an(f)n{h). according to formula [1, (1-5)]. In [1] h is the width of the layer, which we took 
to be h = 1 in our paper without loss of generality. However, the numbers (j)n{h) are not 
known in the inverse problem, since the coefficient n{z) is not known. Therefore formula 
[1, (3.9)] is incorrect. This invalidates the approach in [1]. 

In [1, p. 128] a negative decreasing sequence of real numbers an is defined by equation 
(3.1), which we give for h — 1: 

kVT^n = {n+l)7r + 0{^) (*). 

Such a sequence does not exist: if an < and an has a finite limit then the right-hand 
side of {*) cannot grow to infinity, and if an —00, then the left-hand side of {*) cannot 
stay positive for large n, and therefore cannot be equal to the right-hand side of (*). 
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